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Abstract 

CD , 

^0 ! The properties of wave-like helically twisted normal mode structures on steady 

' relativistic jets are summarized. Wave speeds are a function of the wavelength and 

less than the jet speed. However, normal mode interference can lead to both sta- 
, tionary and superluminal phase effects. A maximum pressure fluctuation criterion 

1/^ ' suggested by numerical simulations of axisymmetric relativistic jets is used to find 

■ the maximum asymmetric jet distortions and velocity fluctuations. Cyclic trans- 

, verse velocity fluctuation can lead to variation in the flow direction on the order of 

\ the relativistic beaming angle. Resulting variation in the Doppler boost factor can 

• lead to significant brightness asymmetries as helical structures twist around the jet 

beam. Growth of these structures is reduced as the jet density, Lorentz factor or 
Mach number are increased. Maximum jet distortion is reduced as the Lorentz fac- 
tor increases and this suggests a reduction in mass entrainment or other non-linear 
disruptive processes that influence the morphological development of radio sources. 
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^ ' 1 Introduction 



Many known or suspected relativistic jets exhibit complex time-dependent 
structures, e.g., 3C 345 (Zensus, Cohen, & Unwin 1995) & M87 (Biretta, 
Zhou, & Owen 1995, Biretta, Sparks, & Macchetto 1999), that are difficult 
to understand in terms of simple static jet bending and/or acceleration of 
components along a fixed trajectory. Wiggling or twisted helical structures on 
multiple length scales provides evidence for wave-like "normal mode" struc- 
tures along these fiows, and precession of the central engine provides a periodic 
driver that can excite the normal modes on multiple length scales. The am- 
plitudes to which normal mode structures can grow has been investigated for 
non-relativistic jets (Hardee, Clarke, & Rosen 1997, Keppens & Toth 1999) 
and for axisjTumetric relativistic jets (Hardee et al. 1998) via comparison be- 
tween theory and numerical simulations. Here we present a prediction of the 
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maximum amplitudes to which asymmetric flow driven (Kelvin-Helmholtz un- 
stable) structures could grow on hydrodynamic relativistic jets and the effect 
on observable jet structure and radio source development. Large amplitude 
normal mode distortions lead to mass entrainment and disruption of highly 
coUimated flow (Rosen et al. 1999, Rosen & Hardee, these proceedings) and 
are an important aspect of the morphological development of radio sources. 



2 The Normal Modes 

A perturbation to a cylindrical jet can be considered to consist of Fourier 
components of the form 

/i(r, 0, z) = /i(r) exp[i{kz incj)- ut)] 

where flow is along the z-axis, and the flow is bounded by r = i? in the 
radial direction. In cylindrical geometry k is the longitudinal wavenumber, n 
is an integer azimuthal wavenumber, for n > the angle of the wavevector 
relative to the flow direction is ^ = tan(n//ci?), and +n and —n refer to 
wave propagation in the clockwise and counterclockwise sense, respectively, 
when viewed outwards along the flow direction. The azimuthal wavenumbers 
n = 0, 1, 2, 3, etc. correspond to pinching, helical, elliptical, triangular, etc. 
normal mode distortions of the jet, respectively. For normal mode n the axial 
wavelength associated with a 360° helical twist of a wavefront around the jet 
beam is given by = nA„ where A„ = 27r//c. 

Propagation and growth or damping of the Fourier components is described by 
a dispersion relation (see Birkinshaw 1991). On the supersonic jet, i.e., Mjt = 
u/ajt » 1 and M^x = u/dex >> 1 where we deflne the sound speeds by = 
FP/ {p+ [r/(r — l)]P/c^} and F is the adiabatic index, each normal mode, 
n, consists of a single "surface" wave and multiple "body" wave solutions to 
the dispersion relation. 

In the low frequency limit all asymmetric (n > 0) growing normal mode 
"surface" wave solutions are given by (e.g. Hardee 1987) 

2 1/2 

u/ku pa — '—^ h ^- 



where rjrei = (aea;/%t)^) and 7 is the Lorentz factor. The growth rate (imagi- 
nary part) decreases and wave speed increases as 7^?7rei >> 1- In the low fre- 
quency limit all "body" wave solutions are given by (see Hardee 1987, Hardee 
et al. 1998) 

kR - k-^R = + 2^ - V2)vr/2 

««~«nm«- ^(M,^2_ 1)1/2 ' 

where n is the mode number, and m > 1 is the integer body mode number. In 
this limit {u!/ku)\reai ~ and the body wave solutions are either purely real 



2 



or damped. At higher wavenumbers, kR > kf^^R > k^^R, the body modes 
are propagating and growing. 

All asymmetric surface and body normal mode waves have a maximum in 
the growth rate (see Gill 1965, Payne & Cohn 1985, Hardee 1987, Hardee et 
al. 1998), at a resonant frequency, wavelength and wave speed: 

uR/aea^ < ul^R/aea, = (n + 2m + l/2)7r/2 , 

27r 7[Mj(2 - m7c2]V2 



[Me^2 _ ^^2/^2] 1/2 + ^[Mj.^2 _ ^^2/^2] 1/2 



U . 



In the above m = indicates the surface wave and m > 1 indicates the 
body waves. At this resonance the growth rate is proportional to {'^Mjf)~^ 
and growth is slowed as ^Mjt increases. 

At high frequencies and short wavelengths the growth rate becomes vanish- 
ingly small and the normal modes propagate like sound waves in the jet with 

l-ajt/u 

u/ku ~ 



1 — ajtu/c^ 



Thus, the normal modes propagate with a wave speed that is a function of 
the frequency relative to the resonant frequency, and wave speeds are always 
less than the underlying flow speed. 



3 Normal Mode Structures 



Displacements of jet fluid can be written in the form 

^(r, 0, z) = A(ro)e*^('-«)e.(i?) exp[t{kzs ± - cut)], 
where Zg and 4>s sue the axial and azimuthal position at the jet surface, Tq is the 
initial radial position, ^r{R) is the radial displacement at the jet surface, and 
the A{r)e^^^'^^ are simple relativistic equivalents (Hardee 1999) of equations 
(AlO) in Hardee, Clarke, & Rosen (1997). The velocity components associated 
with fluid displacements are given by 

Vi(r, (f), z) = d^/dt = —i{uj — ku)^{r, <t), z) . 
The total velocity is v(r, (p, z) = u + vi(r, (p, z) where u is the unperturbed jet 
velocity along the z-axis. The pressure perturbation associated with a fluid 
displacement can be written in the form (Hardee 2000) 

Pi(r, 0, z) = B{rQ)e'^^''°^^r{R) exp[i(A;z, ± tk^s - ut)] . 
A maximal pressure fluctuation criterion, < Fq + -Pi < 2Poi found from 
axisymmetric relativistic jet simulations (Hardee et al. 1998) has been used 
to compute the maximum expected jet distortions for various Lorentz factors 
and Mach numbers (Hardee 2000). Figure 1 shows the maximum normal mode 
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jet cross section distortions for a jet of Lorentz factor 7=10 and with Mjt ^ 
Mf.ri; ^1.3 for the heUcal and elhptical surface wave modes at wavelengths, 
\l'^'*^/R = 8.25, 3.85, 1.88 and A^*'^'*)/^ = 5.54, 3.49, 1.73, where the shortest 
wavelength is at resonance, and for the helical and elliptical first body wave 
modes at the longest unstable wavelength, X^/ R — 3.36 and Af/i? = 2.34. 




-2-1012 -2-1012 -2-1012 -2-1012 




Fig. 1. Maximum cross section distortions for helical and elliptical surface waves 
at three wavelengths (A*, A-', A*), and the first body wave at the longest unstable 
wavelength, A^. 

Jet distortion is a function of the wavelength and scales relative to the reso- 
nant wavelength of the appropriate normal mode and not, for example, to the 
absolute value of the wavelength. High pressure regions are coincident with 
regions where the displacement surfaces are closer together. The cross sections 
rotate counterclockwise as the waves propagate outwards along the jet and, in 
general, high pressure regions precede a trailing low pressure region where the 
maximum outwards surface displacement occurs. For a given pressure fiuctua- 
tion, for example at the resonant wavelength, the accompanying jet distortion 
decreases as the Lorentz factor increases. Thus, in general, the cross sections 
of jets with large Lorentz factors should be less distorted than those with 
small Lorentz factors. Still at long wavelengths the helical and elhptical sur- 
face modes can lead to significant jet distortion. Other results suggest that 
most observable (large scale) structure will develop at wavelengths equal to 
or greater than the resonant wavelengths of the normal modes. Nevertheless, 
tightly wrapped helical structures are possible on hot relativistic jets with 
large Lorentz factors because some resonant wavelengths can be less than the 
jet radius. 

The velocity and pressure fiuctuations that accompany the distortions can 
lead to interesting consequences for Doppler boosted emission features as- 
sociated with asymmetric structures. We can investigate the possible effects 
on emission features by constructing an "apparent" emissivity, P^D"^, around 
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the displacement surfaces shown in Figure 1 where P is the pressure and 
D = {'~f[l — (?;/c) cos^^]}^^ is the Doppler boost factor. Let us assume that 
the initial unperturbed flow angle with respect to the the line of sight is 
^0 = 1/70 = 0.1 rad, the beaming angle. 

Figure 2 gives the values of P'^D^, 7, 9 and D for the surface modes at the 
longest wavelength, A*, shown in Figure 1. Values are plotted around displace- 
ment surfaces with unperturbed radial location of Tq/R = 1/7, 2/7, 7/7. 
The azimuthal angle = lies along the +|/-axis in Figure 1 and increases in 
the clockwise sense. 
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Fig. 2. Emissivity, Lorentz factor, flow angle and Doppler boost factor around the 
displacement contours for the surface wave at the longest wavelength, A*, shown in 
Figure 1. The dashed (dotted) lines correspond to the surface (center) contours in 
the jet cross sections, respectively. 



For the surface waves the highest pressures and highest "emissivities" occur 
at the jet surface and for higher order surface modes the distortion amplitude 
and pressure fluctuation decrease rapidly away from the jet surface. Variation 
between maximum and minimum apparent emissivity as emission enhanced 
regions twist around the jet beam, is about a factor of six for the hehcal mode 
(not shown) and about a factor of three for the higher order elliptical surface 
mode, e.g., the difference in emissivity for the two elliptical filaments in Figure 
2. 

Figure 3 gives the values of P^D^, 7, 9 and D around displacement surfaces 
for the first body modes at the longest unstable wavelength, A^. For the body 
modes the highest pressures and emissivities occur in the jet interior. The 
variations that occur at the resonant wavelength are similar to those shown in 
the Figure 3 for the longest unstable wavelength. Surface mode variations at 
the resonant wavelength are also similar to those shown in Figure 3 although 
the radial structure is different. 
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Fig. 3. Like Figure 2 but for the body wave at the longest unstable wavelength 
shown in Figure 1. 

The asymmetric normal modes lead to high pressure regions helically twisted 
around the jet beam. High pressure regions may be confined close to the jet 
surface, penetrate deeply into the jet interior or be confined to the jet interior 
and range from relatively thin sheets and ribbon like to thicker and tube like 
depending on the mode and wavelength. In general, pressure and emissivity 
maxima are in phase azimuthally in the jet interior at long wavelengths. At 
shorter wavelengths the jet center leads the surface and the high pressure 
region adopts a trailing spiral configuration between jet center and surface. 
Although axial velocity fiuctuation is a small fraction of the relativistic jet 
speed, significant variation in the Lorentz factor occurs on highly relativistic 
jets. Even at high Lorentz factor the velocity fluctuations allow for significant 
angular variation in the flow direction relative to the beaming angle. 



4 Implications 

Jet precession provides a periodic source for the generation of multiple normal 
mode structures. The interference between normal mode structures of differ- 
ent wavelengths and moving at different wave speeds allows for phase effects 
producing superluminal and relatively stationary features. Stationary features 
are created at specific distances from the origin where normal modes of differ- 
ent wavelengths are in phase (Xu, Hardee, & Stone 2000). Superluminal like 
motions have not yet been observed in numerical simulations. Non-relativistic 
numerical simulations have shown resonant pattern motions on the order of 
the resonant wave speed with motions of localized features within about ±15% 
of the average (Hardee 1999, Biretta & Hardee, unpubhshed). 

High Mach number, jet density and Lorentz factor reduce the growth rate 
of normal mode structures and high Lorentz factor reduces the maximum 
amplitude of jet distortions. Nevertheless, significant variation in the Lorentz 
factor and significant angular variation in the fiow direction relative to the 
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beaming angle can occur on relativistic jets. Ultimately, if the jet is sufficiently 
long, large amplitudes at long wavelengths (relative to resonance) can develop 
and lead to mass entrainment, disruption of highly coUimated flow (Rosen et 
al. 1999, Rosen & Hardee, these proceedings), and influence the morphological 
development of radio sources over their hfetime (Eilek, these proceedings). 

The jet in M87 provides an illustration of twisted structures on multiple scales 
that indicate normal mode structures. Radio observations (Owen, Hardee, & 
Cornwell 1989) provide strong evidence for a single bright emission filament 
twisted around the jet beam inside knot D and a fainter twisted filament pair 
between knots D and F. These features indicate the presence of multiple body 
and surface modes. Similar structures are evident in VLBl images at parsec 
scales (Reid et al. 1989). The average motion of well defined structures along 
the M87 jet is subluminal, ^ 0.3c, and should be taken as indicative of a 
wave speed and not the speed of the underlying flow. However, motions less 
than 0.03c and up to about 6c are also observed (Reid et al. 1989, Biretta, 
Zhou, & Owen 1995, Junor & Biretta 1995, Biretta, Sparks, & Macchetto 
1999). Interestingly, the highest observed motions occur for individual bright 
knots often located in regions with complex structure. At HST-1 (Biretta, 
Sparks, & Macchetto 1999) rapidly moving components appear to emerge from 
a subluminally moving component accompanied by brightness variations up to 
±20%. This type of behavior is consistent with interference between different 
wave modes and line of sight phase effects. In the knot A, B and C complex, 
the observed bending and oscillation in jet width may be identified with large 
amplitude helical twisting and elliptical distortion associated with the surface 
modes. This rapidly growing distortion then leads to mass entrainment and 
disruption of the highly coUimated jet flow. 
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